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Abstract 

We introduce a symbolic method for the evaluation of definite integrals containing combinations 
of various functions, including exponentials, logarithm and products of Bessel functions of different 
types. The method we develop is naturally suited for the evaluation of integrals associated with 
specific Feynman diagrams. 
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I. INTRODUCTION 



The elementary integral 
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diverges in the limit v — > — 1, but the singularity can be removed by choosing as arbitrary 
constant 

c : 
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V + 1 ' 

i.e., by subtracting infinities of the same order, in a way to obtain the finite result 



lim L y = In x . 
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This definition is very useful to evaluate integrals containing logarithmic functions. For 
example, its use makes a simple task to prove the following identity (largely used in the next 
sections): 



E„(a,b) = dxx»e~ ax \n(bx) 
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where -0 is the digamma function. For further applications of the definition (jl.2p . see Ref. 



The "trick" of subtracting infinities o: 
of Macdonald function of integer order 



: the same order is at the origin also of the definition 
21 



K n {x) = - hm — , 
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where I u (x) is the hyperbolic Bessel function [3| defined by the series 

(x/2) 2k+v 
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In Ref. [4|, the cylindrical Bessel function has been written as a product of elementary 
functions as follows 
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where the umbral operator c is defined by its action on tp(0) as follows: 

1 



r(/i + i) 
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As a consequence of the relationship 



JJx) = i " iJix) 
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from eq. (II. 6p one has 



x 



I u (x)=^-c) exp jc (^-j | ^(0) 
and, thus, for the Macdonald functions fll.4p the expression 
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As an example, let us consider the case of the function K (x). By applying l'Hopital rule, 
one has 



Kn(x) 
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In (|) /o(x)-exp|c (|) 2 } (lnc)^(O), 



i.e., 



K (x) = - In (|) J (x) - £ (c fc In c) p(0) . 
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Taking into account the definitions (11.21) and (II. 7p . we can write 



(c^ In 6)99(0) = lim 
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and, therefore 



Since 
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(with h = 0) where 7 is Euler's constant, we, finally, obtain 
K (x) = 



j (&rp V2. 



(1.11) 



(1.12) 



(1.13) 



3 



in agreement with the expression reported in Ref. [5|. 

In the following sections the calculation technique outlined in this introduction will be 
used for the evaluation of some definite integrals often encountered in applications, in par- 
ticular in the study of Feynman amplitudes. 



II. INTEGRALS INVOLVING BESSEL FUNCTIONS 



In this section, in order to fix the main results which will be exploited in the second 
part of the paper, we will derive old and new integrals involving Bessel-type functions. The 
starting point will be the general procedure put forward in Ref. {4] to treat the integrals 
of Bessel functions. The technique is fairly efficient and straightforward and benefits from 
the fact that, as already remarked, Bessel functions can formally be treated as elementary 
functions. 

The integral (A > 0) 



dx x^ J\{px) 



-Re A - l<ReyU<-,p>0 



(2.1) 



occurs very often in applications, for example in the evaluation of mass corrections to the 
large square momentum behavior of some QCD processes 6|. Taking into account the 
definition ( II. 6p . we find, in a fairly straightforward way, the well known result (see formula 
6.561.14 in Ref. Q) 
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In the same way, it's easy to show that 
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and therefore, taking into account eq. ( II. 3p . one obtains 
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Moreover, by using the expression (11. lip , we can write (Re fi > — 1) 

M ,o = dxx»K (x) = -jf dxx^ In (|c) exp jc (|) J 



and, by using eqs. (I1.3P and ( ECU]) , 

e*o = 2'- 1 r 2 (£±±) , (2.6) 

that in the case /x = gives 
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The Neumann function (i.e., the function associated with the second solution of the 
cylindrical Bessel's equation) of integer order is defined by a procedure analogous to the one 
leading to Macdonald function, namely 

cos(z/tt) J v (x) - J- V {x) 
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In the case n = one gets 



Y (x) = I Inge) exp|-c (0'} ^(0) . (2.9) 



By using this expression, and eq. (II. 3p . we can easily prove the identity 



vx »Y (x) = ^T 2 (^A Bin(/x|) (-KRe/i<-) 







that, specialized to the case \i = 0, reproduces the well known result 

dxY {x) = 0. (2.11) 
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III. INTEGRALS INVOLVING PRODUCTS OF BESSEL FUNCTIONS 

We have already stressed that integrals of the type (12.11) appears in the solutions to many 
physical problems. In this section we will consider some generalizations of this integral 



associated with specific Feynman diagrams 



0,1s 
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The product of two cylindrical Bessel functions is given by 

j,(-«)4(»«) = (^)" 

with 
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Applying the symbolic method outlined before, this product can formally be rewritten as: 

/ d x\ ^ ( b x\ I/ ( * f X \ ^1 
J^ax) J u (bx) = I— J expl-di^-J \ y^ v {0;a,b) , (3.3) 



where the operator d is defined in such a way that 
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with ip^ v {\\ a, b) generalization to the noninteger case of the function introduced in eq. (13. 21) . 
i.e. 
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By using these results we easily get the following identity (with a and b unequal positive 
number) : 



^a,ii,u(ci'- l b) = / dx x a J^(ax) J y (bx) 
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i.e., 
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This result allows us to evaluate in a very simple way integrals involving products of different 
Bessel-type functions, as, for example, 

/■oo 

~a,n,u( a >fy = / dxx a J fl (ax)Y u (bx) , (3.8) 
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where the definition of Neumann function of noninteger order is obtained removing the limit 
in eq. (I2.8|) . Taking into account eq. (13. 6p . one has 

E a ^ v (a,b) = — j- — - {cos(z/ iv)n a ^ u (a,b) - fi a ft _ v (a, b)} . (3.9) 
sin^ 7i) 

It's evident that the method here outlined can be applied, with a negligible increase 
in calculation complexity, also to integrals involving an arbitrary number of Bessel-type 
functions. 



IV. CONCLUDING REMARKS 



In closing the paper, we apply the operatorial method described in the previous sections 
to the Tricomi-Bessel function 

J„(2Vi) ^ (-!)»*» 

k=0 K ' 

that can be formally written as 

a(x) = c"e-" V (0) (4.2) 

with the operator c satisfying the identity (II. 7p . As a consequence of this equation it's easily 
checked that 



o 
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d*CU*) = -— , / dxC u (x 2 ) = y- (4.3) 



Moreover, from Eqs. (12. 1ft . (12.2ft one has 

> + 1 



^ dss"C„(s 2 ) = V^(2) = l ^ f-KRe//<0. (4.4) 

This function, even though not explicitly mentioned as Tricomi-Bessel function, occurs 
in the evaluation of some Feynman diagrams through integrals of the type 
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where the definition of Macdonald function of noninteger is obtained not taking the limit in 
eq. (jl.4p . By using eqs. flTjjJ) and (13. 6p . we get 

*w(a, 6) = .„ ^ . 7 s [(-1) " fiA^,-,(a, i 6) - ^A-^(a, i &)] . (4.6) 

The variety of topics treated in this paper clearly shows that the method we have proposed 
is particularly efficient to study the problem of the evaluation of integrals of Bessel-type 
functions which very often occurs in physical application. As will be discussed in a forth- 
coming investigation, the intrinsic operational nature of the method make it ideal for an 
automatic implementation by a symbolic manipulator. 
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